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ABSTRACT 


r 


We  consider  the  acceleration  of  solar  flare  protons  by  cyclotron 
damping  of  intense  Alfven  wave  turbulence  in  a magnetic  trap.  The 
energy  diffusion  coefficient  is  computed  for  a near- isotropic  dis 
tribution  of  super-Alfvenic  protons  and  a steady-state  solution  for 
the  particle  spectrum  is  found  for  both  transit-time  and  diffusive 


losses  out  of  the  ends  of  the  trap.  The  acceleration  time  to  a char- 
acteristic energy  i '20  Mev/nucl  can  be  as  short  as  10  sec.  Ow-tH 

_^sis  -ef  phenomenological  arguments  wo^  infer  that  the  Alfven  wave  spec- 
-Z 

trum  has  a u frequency  dependence  and  that  the  correlation  time  of 

'/k-v  '•* 

the  turbulence  lies  in  the  range  5 x 10’S<  < 5 x 10"^  V 


COC^  S ^ cerr  ^ ,Oi  s 


I.  INTRODUCTION 


Solar  flare  protons  are  believed  to  be  accelerated  in  complex 
magnetic  bipolar  configurations  in  active  regions.  Evidence  points  to 
the  idea  that  acceleration  occurs  in  a second  stage  process  [Oe  Jager, 
1969]  in  contrast  to  primary  energy  release  of  the  flare  in  a first 
stage  process  which  is  manifest  directly  in  quasi-thermal  heating  and 
non-thenral  10-500  Kev  electrons  [Svestka,  1976].  The  association  of 
Type  II  radio  bursts  with  proton  flares  is  compatible  with  the  scenario 
that  a noisy  MHO  shock  is  generated  soon  {^3  mins)  after  the  initiation 
of  the  flash  phase  and  accelerates  protons  by  stochastic  processes  as  it 
propagates  upwards  and  away  from  the  flare  site  [McLean  et  al.,  1971; 

Lin  and  Hudson,  1976].  The  shock  is  certainly  energetic  enough  and 
current  belief  is  that  Fermi  processes  are  at  work  in  conduction  with 
wave-particle  scattering  to  energize  protons  up  to  energies  > 10  Gev 
[Melrose,  1974;  Svestka,  1976].  A large  part  of  our  theoretical  under- 
standing of  acceleration  mechanisms  has  derived  from  efforts  directed 
towards  understanding  the  energization  of  galactic  cosmic  rays  [Rosen, 
1965]. 

In  this  paper  we  will  attempt  to  model  the  solar  proton  acceleration 
process  with  a very  simple  field-plasma  configuration  which  might  be 
described  as  a "leaky  trap".  Energization  takes  place  in  the  corona 
in  closed  magnetic  field  structures  but  allowance  is  made  for  particle 
losses.  As  a first  step  in  this  direction  we  have  considered  only 
particle  escape  along  the  flux  tube  to  the  low  solar  atmosphere.  The 
acceleration  mechanism  is  cyclotron  damping  by  super-Alfvenic  protons 


of  intense  Alfven  waves  propagating  in  the  trap.  This  process  is 
intrinsically  stochastic  and  no  adiabatic,  Fermi-type  mechanism  is 
invoked  for  the  acceleration. 

Davis  [1956]  gave  the  earliest  description  of  the  effect  that  a 
turbulent  state  of  MHD  waves  produces  a diffusion  in  pitch  angles 
towards  isotropy.  The  development  of  quasi -linear  theory  [see  Kennel 
and  Englemann,  1966,  and  their  refs.]  provided  a useful  analytic  tool 
with  which  one  could  address  the  problem  of  diffusion  in  phase-space 
resulting  from  turbulent  wave-particle  interactions.  Sturrock  [1966], 
using  a Fokker-Planck  formalism,  achieved  similar  results  and  thus 
fortified  our  confidence  in  these  methods. 

Melrose  [1974]  demonstrated  that  the  combined  effects  of  stochastic 
pitch  angle  scattering  by  Alfven  waves  and  Fermi  acceleration  could 
account  for  the  energization  of  solar  flare  protons.  The  source  of  the 
Alfven  turbulence  was  considered  to  be  the  non-thermal  protons  them- 
selves which,  if  sufficiently  anisotropic  in  pitch  angle  distribution, 
would  drive  an  electromagnetic  Alfven  wave  instability.  Since  the  free 
energy  for  the  waves  is  derived  from  the  kinetic  energy  of  the  non- 
thermal  protons  in  this  picture,  another,  independent,  free  energy 
source  must  energize  the  protons  (e.g.,  a Fermi  process).  In  our  model 
the  energy  flow  is  opposite  to  that  of  Melrose.  We  assume  that  Alfven 
turbulence  is  generated  ad  hoc  by  an  exciting  agency  that  is  flare 
related.  Using  this  approach  we  may  eliminate  the  Fermi  mechanism  and 
consider  the  acceleration  possible  from  just  short  wavelength,  kv  = n^p. 
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Alfven  waves  which  not  only  scatter  in  pitch  angle  but  also  energize 
through  collisionless  damping  [Schatzman,  1967;  Sturrock,  1974]. 

The  model  can  be  applied  to  the  system  described  by,  for  instance, 
McLean  et  al.  [1971]  in  which  a localized  shock  generates  Alfven  waves 
propagating  throughout  a flux  tube.  With  some  modifications,  the  model 
may  also  be  applied  to  a system  like  that  of  Sturrock  [1974]  in  which 
protons  are  accelerated  on  open  field  lines  by  a noisy  MHD  shock. 
However,  recent  Y-'"ay  observations  have  suggested  that  some  protons  may 
be  accelerated  as  early  as  the  impulsive  phase.  We  have,  therefore, 
concentrated  on  the  consequences  of  a leaky  trap  without  pinning  down 
the  source  of  the  turbulence  or  when  and  how  it  is  generated.  If 
impulsive  phase  proton  acceleration  is  confirmed,  then  we  hope  our 
model  may  provide  insight  into  the  primary  energy  release  mechanism 
itself. 

In  Section  II  we  outline  the  ingredients  of  the  model  and  inquire 
into  the  possibility  of  proton  acceleration  during  primary  energy 
release.  In  Section  III  we  derive  the  energy  diffusion  coefficient, 
^EE*  cyclotron  damping  mechanism  for  a near- isotropic  distribu- 

tion of  super-Alfvenic  ions.  In  Section  IV  we  apply  this  result  to 
the  leaky  trap  with  either  ballistic  transit-time  losses  or  diffusive 
losses  out  the  ends  of  the  trap.  In  the  final  section  we  summarize 
our  results  and  comment  on  the  possibility  of  flare-generated  MHD 
standing  waves  in  the  trap. 
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II.  MODEL 

A.  Basic  Assumptions 

The  basic  magnetic  field  configuration  we  shall  consider  Is  that  of  a 
closed  loop,  bipolar  arch  structure.  We  presume  that  some  exciting 
agency,  that  Is  flare  related,  generates  Intense  MHD  turbulence  which 
Interacts  with  particles  at  coronal  heights.  The  closed  field  lines 
serve  as  a magnetic  trap  for  the  non-thermal  particles  which  are  stochas- 
tically accelerated  by  the  turbulence  but  are  prevented  from  escaping  by 
mirror  forces  arising  from  the  convergence  of  the  mean  field  towards 
the  photosphere.  The  turbulence  Is  assumed  to  consist  of  short  wavelength 
perturbations,  ^ « L,  much  smaller  than  any  macroscopic  scale  of  the 
system;  this  enables  us  to  use  homogeneous  quas1-l Inear  theory  to  gauge 
the  acceleration  efficiency.  The  dominant  Interaction  with  the  waves  Is 
considered  to  be  with  the  wavevector  component  parallel  to  [I.e., 

o]  and  we,  therefore,  model  the  turbulence  as  an  ensemble  of  parallel- 
propagating  Alfven  waves  In  the  frequency  range  Wg  < w « a^p,  where  n^p 
Is  the  proton  gyrofrequency  and  Is  a low  frequency  cutoff. 

The  question  of  particle  phase-space  losses  Is  an  Important  one  and 
different  assumptions  concerning  what  loss  mechanism  prevails  will  alter 
the  picture  significantly.  We  shall  consider  proton  energies  sufficiently 
large  that  they  may  be  considered  collisionless  with  respect  to  binary 
coulomb  Interactions  [Braginskil,  1965]  during  the  acceleration  period. 
However,  If  the  protons  remain  In  dense  regions  sufficiently  long,  this 
deenergization  process  must  be  accounted  for.  We  may  expect  that  some 
cross-field  diffusion  of  particles  out  of  the  trap  Is  likely  to  occur 
during  the  acceleration  period.  However,  we  have  chosen  to  consider  what 
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may  be  the  fastest  loss  process:  particle  escape  along  the  field  to  the 
chromosphere  and  photosphere.  A consequence  of  our  assumption  of  intense 
Alfven  turbulence  (with  a predominant  magnetic  component)  is  that  particles 
will  be  scattered  very  efficiently  in  pitch  angle  and  can  either  stream 
out  or  diffuse  out  of  the  acceleration  region  along  the  flux  tube  to  the 
lower  atmosphere.  For  either  case  we  have  scaled  the  losstime  to  the 
particle  ballistic  travel  time  and  the  diffusion  time  along  the  length  of 
the  flux  tube.  In  the  first  case  (transit-time  losses)  our  model  is  anal- 
ogous to  the  process  described  by  Kennel  [1969]  who  considered  stochastic 
acceleration  of  electrons  in  the  earth's  magnetosphere.  In  the  second 
case  (diffusive  losses)  the  model  incorporates  diffusive  spatial  transport 
in  the  presence  of  elastic  wave  collisions  [Jokipii,  1971]  analogous  to 
the  propagation  effect  observed  in  interplanetary  space. 

B.  Applications 

The  association  of  Type  II  radio  bursts  with  proton  events  observed 
at  earth  suggests  that  a flare-initiated  shock  propagates  upwards  and 
accelerates  protons.  If  the  shock  is  the  exciting  agent  of  MHD  waves, 
then  our  model  may  be  used  to  compute  the  acceleration  of  protons  in  a 
configuration  like  that  proposed  by  McLean  et  al.  [1971].  It  must  be 
determined,  for  instance,  if  the  acceleration  takes  place  ahead  of  or 
behind  the  shock,  and  the  question  of  how  the  protons  obtain  access  to 
interplanetary  field  lines  must  be  addressed.  If  the  field  lines  of  the 
trap  "open  up"  [Kopp  and  Pneuman,  1976], then  particles  will  be  injected 
ball istically  into  the  solar  wind;  but  if  the  field  lines  remain  closed. 


but  distended,  cross-field  diffusion  effects  must  occur  to  exhaust  the 


trap.  Chronologically,  this  process  is  an  "after  effect"  of  the  flare; 
that  is,  subsequent  to  or  during  the  primary  energy  release  the  shock 
is  formed  and  proton  acceleration  begins  at  times  ^ 3 mins  after  the 
initiation  of  the  flash  phase. 

However,  y-ray  observations  [Chupp  et  al.,  1973]  made  by  the  OSO-7 
satellite  have  indicated  that  proton  acceleration  may  begin  at  times  coin- 
cident with  the  flash  phase  onset.  Figure  1 is  a plot  of  the  Aug.  4,  1972 
2.2  Mev  Y-ray  line  together  with  hard  x-ray  profiles  from  Van  Beck  et  al. 

[1973].  We  note  that  while  the  maximum  of  the  y-ray  flux  occurs  after 
the  hard  x-ray  maximum  a not  insignificant  flux  is  measured  very  soon  after 
onset  [Svestka,  1976].  In  fact,  the  flux  at  0625  UT  is  1/3  of  the  maximum. 

Since  there  must  be  a finite  acceleration  time  for  the  protons  to  reach 
these  energies  in  sufficient  number.  Figure  1 indicates  that,  within  the 
temporal  resolution  of  the  instrument  and  ignoring  any  nuclear  reaction 
time  delay,  a generous  upper  limit  for  the  acceleration  time  -v-  3 mins 
(assuming  that  energization  begins  at  0622).  The  acceleration  time  could 
be  much  less  (10s  < < 200s)  and  still  be  consistent  with  the  non- 

thermal  2.2  Mev  emission  measure  reaching  a maximum  later  as  more  flux 
tubes  become  involved  and  the  acv-releration  region  increases.  The  measure- 
ment at  0622,  if  real,  gives  support  for  a short  '\'10s  acceleration  time. 

Ramaty  et  al.  [1975]  interpret  a break  in  the  hard  x-ray  spectrum  of 
Aug.  4,  1972  at  about  500  Kev  towards  a harder  profile  at  higher  energies 
as  evidence  for  electrons  undergoing  second  stage  acceleration.  t 


Since  we  might  expect  intense  Alfven  wave  turbulence  generated  in 
primary  energy  release  mechanisms  like,  for  instance,  magnetic  field 
tearing  instabilities  [Spicer,  1976],  we  have,  therefore,  considered 
a parameter  regime  for  the  leaky  trap  that  might  be  applied  to  plasma- 
field  conditions  in  the  vicinity  of  x-ray  kernals  and  their  coronal 
loops  at  the  time  of  the  impulsive  phase.  The  model,  then,  can  allow 
us  to  infer  something  about  the  primary  energy  release  process  itself 
rather  than  post-explosion  conditions  in  MHD  shocks.  Future  observations 
should  define  clearly  the  role  of  magnetic  turbulence,  especially  with 
regard  to  primary  energy  release,  and  can  test  our  supposition  that 
proton  acceleration  can  occur  as  early  as  the  impulsive  phase. 


If 


III.  DIFFUSION  COEFFICIENT 
A,  Quasi -I inear  Theory 

We  consider  a homogeneous  plasma  in  a magnetic  field  directed 
along  the  z axis  and  assume  that  the  turbulence,  composed  of  parallel- 
propagating  Alfven  waves  (k^  = o)  with  frequencies  in  the  range 

o^o  < 0)  = |k|V^  « (1) 

can  be  described  by  a spectral  density 

< ^I/Ott  > = 2 /“du.4(o)).  (2) 

O 

Va  = BQ/(4TrMnjj)**  is  the  Alfven  velocity,  n^p  = eB^Mc  is  the  proton 
gyrofrequency,  and  is  the  perturbation  magnetic  field  related  to  the 
induced  electric  field  by 

(3) 


The  factor  of  two  in  equation  (2)  is  to  discriminate  between  waves  travel- 
ling parallel  and  anti  parallel  to 

It  is  well  known  that  the  temporal  evolution  of  a distribution  of 
particles,  say  protons,  in  response  to  resonant  scattering  by  Alfven  waves, 
can  be  described  by  a quasi-1 inear  relation  [Kennel  and  Engleman,  1966; 

Hall  and  Sturrock,  1967;  Davidson,  1972]  which,  for  kx  = o,  reduces  to 
in  the  non- relativistic  regime 

3 J 3 


where 


Air^e^Vl  |Vg|£(a.) 
M2c2  IVg-v^l 


+ 
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In  the  frequency  range  consideredf  the  group  velocity  is  the  Alfven  speed, 
|Vg|  = V^;  and  are  unit  vectors  in  velocity  space  (cylindrical 
coordinates),  azimuthal  symmetry  about  being  assumed.  The  effect  of 
the  finite  perturbation  electric  field  (3)  is  represented  in  the  terms 
consisting  of  the  integer  1.  The  unique  relationship  between  the  proton 
and  the  Alfven  wave  in  this  interaction  is  expressed  in  the  cyclotron 
resonance  condition  which  is  given  as 

“ - '"''z  = -“cp  ’ 

the  upper  (lower)  sign  taken  for  the  L (R)  mode  of  the  Alfven  wave.  For 
k > 0 the  terms  L and  R also  give  the  polarization  of  the  wave.  Equation 
(6)  indicates  that  an  effective  "collision"  occurs  when  the  doppler- 
shifted  wave  frequency  is  the  proton's  gyrofrequency  and  the  sense  of 
rotation  is  the  same. 

For  the  range  of  parameters  we  shall  consider,  « cand  the  effect 
of  the  predominantly  magnetic  wave  is  to  scatter  the  particles  primarily 
in  pitch  angle.  It  is  convenient,  therefore,  to  express  (5)  in  spherical 
coordinates.  If  cose  = v^/v  we  have,  using  (6), 


M2c2 


V..-  V, 


kVx  1 r - 
sine  ± TT — cose  I — + cose  + ^ — si 
cp  cp 


ne 


1 3f  1 r ^ r - ) - 

^ ^ I sine  ± — cose  v + cose  + ^ — sine  e . 
''  l^i  «cp  i i “cp  ^ 


(7) 


Inspection  of  (6)  indicates  that  if 


01  « then  | v^|  » if 


resonant  scattering  is  to  occur.  Protons  with  |v 


< follow  the  wave 
mv^ 


adiabatically  conserving  their  first  invariant  y = during  the 
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oscillation,  whereas  super-Alfvenic  protons  can  satisfy  (6)  and  are 
non-adiabatically  scattered,  the  interaction  being  sufficiently  rapid 
(t  n~p)  so  as  to  violate  their  first  invariant.  Equation  (6)  then 
permits  a useful  perturbation  expansion  [Melrose,  1974] 

+1/VZ  ■ I 

l-oi/kv. 


k/fi 


cp 


= 1 uu  I uu  u I VII  cApnii^  I VII  ^ric  i i ' 

■ ] 


(8) 


for  particles  with  velocities  |v^|  » V^.  Retention  of  the  zero  order 
term  in  (8)  provides  strictly  magnetic  effects,  while  the  first  order 
term  incorporates  effects  of  a finite  electric  field  due  to  ui  f o.  The 
expansion  of  (7)  to  zero  order  yields 


r(o) 


47r2e2 


1 


3f  I 
30  ® 


(9) 


M^C^  |COS0 

Any  reference  to  frequency  here  is  incidental  since  if  (2)  is  given 
alternately  as 

< > =J'“dkE(k)  (10) 

we  have  the  relation 

V^fi(a))  =fi(k)  (11) 

and  thus  (9)  can  also  represent  non-adiabatic  pitch  angle  scattering  from 
stationary  magnetic  inhomogeneities  described  byfi(k)  [Jokipii,  1971]. 

On  the  shortest  time  scale,  the  turbulence  acts  to  isotropize  the 
distribution  and  on  a longer  time  scale  0(v^/vp  the  particles  are  ener- 
gized through  stochastic  accelerations  from  This  results  in  a diffu- 
(2) 

Sion  current  J*  ' which  we  may  conveniently  write  using  the  approximation 
that  the  distribution  is  nearly  isotropic  from  pitch  angle  scattering  and 
^ = 0.  Retaining  the  first  order  term  in  (8)  provides  the  second  order 
term  for  ^(2)  ^ 0(VS/v2) 
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A 


r 


3(2) . 

v|cos0| 

and  the  distribution  evolves  as 


v|l 

M^c^-  loose  I v2  3v  L 


(12) 


(13) 


B.  Wave  Spectrum 

To  inodel  the  turbulence  we  assume  the  spectral  density  can  be 
represented  by  a power  law  in  frequency 

S(o))  - Bq  ^ ^ 0(<^Wq)  . (14) 

9(x)  is  the  Heaviside  step  function  which  introduces  a low  frequency 
cutoff  at  Up.  It  is  also  implicit  in  (14)  that  frequencies  do  not  extend 
as  high  as  the  proton  gyrofrequency.  The  resonance  condition  (6)  gives 
the  inverse  wave-particle  relation 


U)  = |k|V^  = P. 


(15) 


vicosel 

The  subsequent  steps  are  straightforward:  the  spectral  density  is 
normalized  to  the  perturbation  field  energy  density  in  (2)  to  determine 
equations  (14)  and  (15)  are  inserted  into  (13)  and  the  equation  is 
averaged  over  pitch  angles,  and  the  change  of  variable  to  energy  per  nucl 
is  made. 

Generalizing  this  procedure  for  an  arbitrary  ion  of  ionic  charge  Ze 
and  mass  AM,  the  temporal  evolution  of  the  differential  number  density  for 
any  ionic  species  can  be  written  as 


(16) 
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where  the  total  number  density  of  the  ionic  species  is  given  as 


N = /“dE  N(E)  (17) 

and  D = (A/Z)^""^  Dp  , (18) 

B?/87r  2n 

= Znn  -! ^ ^ (19) 

P ‘^P  B2/87t  (2n+l)(2n+3)  Ejjp 

2 2 ®o 

The  variable  E = ^sMv^  = energy/nucl  of  the  species,  E^  = 

is  the  Alfven  energy,  M = proton  mass,  and  the  index  n is  related 

to  the  wave  spectral  index  by 

n = i4(6-l)  . (20) 

The  normalization  of  A(a))  introduced  the  low  frequency  cutoff 

into  the  scaling  of  the  diffusion  coefficient  and  we  have  written  it  in 


energy  units  as 


^mp  = 


( ^ ]‘ 
I 


Emp  is  related  to  the  maximum  energy  a proton  can  achieve  in  the  following 
way:  if  E^^^p  is  non-relativistic  in  value,  then  p = E^j^p,  but  if  E^^^p 
is  relativistic,  one  must  use  the  relativistic  form  of  the  resonance 
condition  instead  of  (6)  to  determine  E„.„  Thus,  the  maximum  energy 

inaX  f p 

a proton  can  be  accelerated  to  is 

^ : E„_  « IGev 

I mp  mp 

^max.p  ' ''a  "cp 

f » tGPV 


Again,  so  long  as  we  are  considering  the  non-relativistic  regime  of 
energies  in  (16),  (21)  is  used  by  definition  for  all  values  of  “q. 

It  is  of  interest  to  ask  what  minimum  energy  (injection  energy)  of 
electrons  is  necessary  for  them  to  be  accelerated  by  Alfven  waves.  The 
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resonance  condition  for  electrons  with  gyrofrequency  n = eBp/mc  is 


given  as 


which  indicates  that  if  w < we  have  [Sturrock,  1974] 

Elnj.p=  Ea  = 


^inj.e 


non-rel 


Mc^  — : ultra-rel 


These  restrictions  apply  for  MHD  Alfven  waves  where  Extension 

of  the  wave  spectrum  to  higher  frequencies  to  include  whistler  waves 
(R  mode)  or  Langmuir  waves  lowers  the  electron  energy  threshold  [Melrose, 
1974],  but  is  better  addressed  as  a separate  problem  since  the  wave 
properties  are  quite  different  and,  phenomenologically,  acceleration  by 
higher  frequency  waves  may  be  related  to  first  stage  acceleration 
processes  [Hoyng,  1977]. 

If,  instead  of  energy,  one  considers  a diffusion  in  the  non-relativis- 
tic  regime  of  rigidity,  R = AMvc/Ze,  we  have  alternatively 

3 3 r n2n+2  3 


^ N(R,t)  = K ^ ^ (N/R2)  ] 


where 


1/  — A IT  _ 

■ I 7 “cp 


B2/0.  (2n+l)(2n+3)  R' 


C.  Physical  Interpretation 

An  interpretation  of  the  physical  process  leading  to  equations  (4), 
(6),  and  (14)  is  given  in  Figure  2.  We  have  decomposed  the  turbulence 
into  an  ensemble  of  elemental  scattering  centers  of  length  x < 

This  ensemble  is  equivalent  to  the  actual  situation  in  that  the  idealized 
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scatterers  have  the  same  autocorrelation  function  as  the  actual  system, 

p(t)  ] (27) 

p(x)  -.i"dk|(k)e^''’‘/[i“dkA(k)  ] . (28) 

From  the  power  spectrum  given  in  (14)  we  identify  (6  > 1) 

corr  corr  A corr  ' ' 

and^(l()  [or£(  (d)]  is  a measure  of  the  differential  number  density  of 

elemental  scattering  centers  of  size  k'^  -v  x < x . A super-Alfvenic 

corr 

proton  satisfying  the  resonance  condition  (6)  can  be  effectively  scattered 
as  drawn  in  Figure  2.  Particle  A sees  a component  of  parallel  to  v 
throughout  the  interaction  and  is  accelerated,  while  particle  B is  180® 
out  of  phase  and  is  decelerated.  A collection  of  random  (cyclotron) 
phased  particles  will  experience,  however,  a dispersion  in  energies. 

The  acceleration  process  described  here  cannot  be  interpreted  as  a 
stochastic  Fermic  type.  The  acceleration  mechanism  described  by  Fermi 
[1949,  1954]  involves  an  adiabatic  process  such  that  the  particles'  first 
and  second  invariants,  y = '’'''V2B,  J =y*V||ds,  are  conserved.  Pitch 
angle  scattering  is  an  ancillary  process  invoked  to  counterbalance  the 
systematic  decrease  in  the  particle's  pitch  angle  as  it  is  accelerated 
[Davis, 1956].  We  have  described  a process  where  oscillations  are  suffi- 
ciently rapid  that  they  non-adiabatically  violate  the  first  invariant. 
Energization  (or  denergization)  is  concomitant  with  pitch  angle  scattering 
when  0)  0 and  such  a process  is  best  interpreted  as  cyclotron  damping 

of  Alfven  waves  [Stix,  1962]. 
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IV.  CHARACTERISTIC  SCALINGS 


In  order  to  apply  the  result  of  equation  (16)  to  the  leaky  trap 
model,  we  have  considered  three  problems:  a transient  solution  and  a 
steady-state  solution  with  transit-time  and  diffusive  losses. 

A.  Initial  Value  Problem 

Two  facts  which  immediately  emerge  from  injection  of  (16)  is  that 
the  energization  time  scales  as  -v-  E^"^  and  that  the  solution  asymptotes 
to  a "steady  state"  N(E)  E"*'.  The  energization  time  indicates  that 
n = 1 distinguishes  two,  quite  different,  diffusion  regimes  which  favor 
acceleration  of  either  high  or  low  energy  particles.  We  shall  argue  in 
this  paper  that  the  regime  which  favors  low  energy  particle  acceleration 
is  more  appropriate  for  solar  proton  flares  and  we  shall  consider  the 
parameter  regime 

1<6<3  ■*-*■  0<n<l.  (30) 

This  statement  is  equivalent  to  the  expectation  that  if  the  power  spec- 
trum in  (14)  is  sufficiently  flat,  protons  can  be  accelerated  over  a 
large  energy  range  in  the  shortest  time  for  a given  perturbation  field 

i 

energy  density.  Increasing  6 to  larger  values  puts  all  the  energy  in 
the  low  frequencies  and  biases  the  acceleration  in  favor  of  very  high 
energy  particles  [Equation  (15)]. 

The  homogeneous  solution  of  (16)  with  the  boundary  conditions  at  an 
injection  energy  E^  that 

N(EQ,t)  = Nq  , N(E>EQ,t=0)  = 0,  N(E  -^  -)  -^  o 


j 

i 

i 

A 

L 


r 


is  given  as 
N(E 


J” 

t)  - 

-ioo 


e'  ( I ^ 


kJee"  A') 


with 


P = »t(l-2n) 
w = Ss(l-n)  > 0 


^0  K (ceJJ  /s) 
u 0 ' 


_ 2n+l 

■ ?rr??y  • 


(31) 


(32) 


If  we  consider  times  long  enough  such  that^E^  « 1,  the  small- 
argument  form  of  K (z)  can  be  used  in  the  denominator  and  the  integral 


is  evaluated  as 


N(E.t)  ■ N„  [-^ 


E„  r(y,  HC^E^w) 


r(u) 


(33) 


where  r(a,x)  is  the  incomplete  gamma  function  [Abromowitz  and  Stegun,  1970]. 
This  function  exhibits  two  different  behaviors  depending  on  whether 


X = 2 e2W  > 1 _ 

We  may  approximate  the  solution  as 


(34) 


N, 


Eq  .n 

T 


X <<  1 


N(E,t)  = 


X » 1 


r(„) 

The  solution  changes  from  a power  low  at  low  energies  to  an  exponential 
behavior  (w>0)  at  the  diffusion  front  characterized  by  a break  energy. 
Eg,  which  depends  on  time  and  is  defined  as  x = 1 leading  to 


(35) 


1 

EB(t)  = [(l-n)2  Dt]  ^ 


(36) 


r 

I 


In  a subsequent  section  we  shall  argue  for  the  feasibility  of  the 
value  n = Anticipating  this  discussion  and  defining  £g  = 
scale  the  break  energy  as 


I A,  r I i )■  I le  1 1 ) 


lOGev 

"^p 


t^  Mev/nucl  . 

Time  is  in  units  of  seconds  and  our  choice  of  and  n. 


(37) 


.Q  — ..Q  background 

density  corresponds  to  an  Alfven  energy  of  ■ 41.4  Kev.  Note  also 

[Equation  18]  that  for  n * *j.  reference  to  7/A  disappears  and  (37)  is  valid 

for  any  ionic  species.  The  value  10  Gev  for  E corresponds  to  E , » 

mp  max  p 

4.3  Gev. 


B.  Steady- State  Solution 

Assuming  that  proton  acceleration  can  occur  for  a period  longer  than 
10  sec,  (37)  admits  the  possibility  of  very  efficient  proton  acceleration. 
Eg  being  a measure  of  the  "hardness"  of  the  proton  spectrum.  A realistic 
model  should  allow  for  particle  losses  and  phenomenologically  we  may  write 


(38) 


to  describe  the  time  asymptotic  evolution  of  the  transient  solution  (33). 
The  steady-state  solution  of  (38)  should  give  a meaningful  result  provided 
that  the  fastest  losstime 


T 


(39) 
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can  be  ascertained.  The  solution  for  q = Js(l-n-e)  > 0 is 


N(E)  = N ^ , 

° ^ ^0  J K (XE^) 

V 0 


where 


P = J6(l-2n) 


TTiHrel 


i5r 


q = Js(l-n-B)  > 0 


If  xeJ  <<  1 we  may  use  small  argument  form  for  the  denominator  and 


(40)  exhibits  the  following  behavior 


N(E)  = 


No  (Eo/E)^ 


: XE*’  « 1 


N f JsXE 

° r(v)  ' 


J ^ Jit(n.e)  o-XE'’  . ^pq^^i  ^ 


The  steady-state  break  energy,  defined  as  xE^  = 1,  is 


Eg  = [ Dto  ] . 


C.  Transit-Time  Losses 

In  a magnetic  field  minor  geometry,  particles  with  pitch  angles 

less  than  the  loss  cone  value 

sin^e  = B /B  (44) 

"o  o'  max  ' ' 

can  penetrate  to  a field  strength  of  where  they  may  be  considered  as 

lost  [through  significant  denergization  processes  in  the  dense  chromosphere 

and  photosphere].  This  angle  corresponds  to  a loss  cone  "volume"  of 

4TT[l-coseo]  steradians.  If  pitch  angle  scattering  occurs  rapidly  enough 

to  maintain  a significant  number  of  particles  in  the  loss  cone,  then  the 
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losstime  is  estimated  to  be  the  transit  time  along  the  flux  tube  with  a 
weighting  factor  of  the  fractional  volume  of  the  loss  cone  [the  strong 
pitch  angle  diffusion  regime  of  Kennel,  1969],  i.e., 

T = {n/Z)^  L(l-cosejj)‘^  E"**  . (45) 

Transit-time  losses  lead  to  a value  of  3 = -*s  in  (39).  We  note  that 
if  n = is,  then  the  solution  (40)  behaves  at  large  energies  as  N(n=-3=is) 

E ^ e . Observations  made  at  lAU  have  indicated  that  over  a large 
range  of  energies,  proton  spectra  can  be  fit  by  an  exponential  in  rigidity, 
the  fit  improving  at  high  energies  (^200MV  = 21Mev)  [Freier  and  Webber, 

1963;  Van  Hollebeke  et  al.,  1975].  An  exponential  in  rigidity  is  also 
consistent  with  models  based  on  y-ray  observations  [Ramaty  and  Lingenfelter, 
1972].  On  this  basis  we  argue  for  the  feasibility  of  the  value  n = *$ 
corresponding  to  a wave  spectral  index  of  6 = 2. 

If  we  take  a coronal  field  strength  of  1006  and  a photospheric  field 
of  10006,  the  loss  cone  angle  = 18°.  Using  (45)  in  (43)  gives  a steady- 


state  break  energy  for  transit-time  losses  of 


Although  (46)  gives  a reasonable  value  for  the  turnover,  the  scaling 
may  not  be  correct.  In  the  absence  of  cross-field  losses,  (45)  without 
the  loss  cone  weighting  factor  must  be  the  shortest  conceivable  time  to 
deplete  the  trap.  But,  for  the  parameters  chosen,  the  mean  free  path  for 


19 


pitch  angle  scattering  is  much  less  than  the  length  of  the  tube.  Under 
these  circumstances,  a spatial  diffusion  along  the  field  to  the  lower 
atmosphere  may  be  more  appropriate  than  losses  scaled  to  the  particle 
ballistic  transit  time. 


D.  Diffusive  Losses 

If  the  phase-space  distribution  function  varies  along  the  z direction 
and  we  consider  only  the  effects  of  pitch  angle  scattering  by  waves,  then 
(4)  in  zero  order  can  be  written  as 


|i  * ^ f (z.f.t)  = ^ ] (47) 

and  n = cose.  Jokipii  [1971]  has  shown  that  in  the  approximation  that  f 
is  nearly  isotropic,  (47)  reduces  to  a diffusion  equation  for  the  differ- 
ential number  density 

^N(z,E,t)  = K,,  ^ , (48) 

at  II  ^^2 


where 


■'119 

and  from  (9)  we  have  for  protons 


i C/'X''  “nn  ]■’ 


""  „2c2*^'  ^ V 


nn 


(49) 


(50) 


Substituting  the  spectral  density  (14)  in  the  above  and  generalizing  for 
an  arbitrary  ion  yields  the  result 


/ 

-1 


1 


dn  0 


nn 


^1/8n 

B2/8ir 
0 


2n  e" 
(^n+l)(2n+3)  nn 

^mp 


(51) 
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and  the  parallel  spatial  diffusion  coefficient  is 


V.  SUMMARY  AND  DISCUSSION 

We  have  considered  some  of  the  consequences  of  a turbulent  spectrum 
of  intense  Alfven  waves  propagating  along  the  magnetic  field  of  coronal 
flux  tubes.  The  model  describes  a steady-state  situation  where  a popula- 
tion of  ions  is  maintained  at  a threshold  velocity  of  the  Alfven  speed 
and  is  scattered  in  pitch  angles  by  the  turbulence  towards  isotropy 
and  energized  concomitantly  through  collisionless  damping.  The  accelera- 
tion time  to  -vlOMev  energies  can  be  as  short  as  'v-lOs  and  the  process  itself 
is  limited  only  by  how  well  the  trap  can  contain  the  accelerated  particles. 
The  measure  of  the  acceleration  efficiency  has  been  gauged  by  a break 
energy  of  the  particle  spectrum  which  characterizes  a turnover  from  a 
relatively  flat  partical  spectrum  to  a rapidly  decreasing  exponential  at 
. high  energies.  This  break  energy  is  approximately  the  value  at  which  the 

acceleration  time  equals  the  loss  time. 

The  primary  loss  mechanisms  which  we  have  considered  are  transit- 
time losses  and  diffusive  losses  out  of  the  ends  of  the  flux  tube.  For 
the  range  of  parameters  such  that  the  mean- free- path  for  pitch  angle 
scattering  is  small  relative  to  the  length  of  the  tube,  the  diffusive 
picture  is  probably  more  correct.  The  question  of  other  losses  is  impor- 
tant and  further  work  is  needed  to  clarify  the  role  of  cross-field  diffu- 
sive effects,  especially  with  regard  to  how  particles  obtain  access  to 
interplanetary  field  lines.  However,  if  the  acceleration  region  consists 
of  open  field  lines  [Sturrock,  1974]  this  problem  disappears. 

I . • One  important  conclusion  we  have  drawn  is  that  for  a given  r.m.s. 

I magnetic  perturbation  B-j  the  wave  power  spectrum  in  (14)  must  be  relatively 
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flat  (30)  if  efficient  acceleration  is  to  occur.  For  large  values  of 
the  wave  spectral  index  6 the  diffusion  process  favors  the  acceleration 
of  the  very  high  energy  particles  which  are  scattered  by  the  most  ener- 
getic oscillations  near  the  correlation  length  of  the  turbulence;  in 
this  case  the  energy  diffusion  coefficient  scales  as  a strong  power  of 
(Eo/Emp)'^  [Newman,  1975]  indicating  that  the  acceleration  of  low  energy 
particles  near  the  injection  energy  is  inefficient  and  holding  up 
the  works.  On  the  basis  of  this  model  proton  flares  are  associated  with 
both  large  magnetic  field  perturbations  and  flat  power  spectra  (30). 

This  result  is  consistent  with  the  findings  of  La  Combe  [1977]  who  con- 
sidered particle  acceleration  in  radio  galaxies  by  Alfven  waves.  We  may 
also  extend  the  result  (19)  to  include  the  value  6 = l(n=0)  for  which 
the  substitution  is  made  in  both  (19)  and  (51) 


6 = 1 3 


ln“max/“o 


In  particular,  we  have  shown  that  the  value  6=2  leads  to  an  expo- 

-k  -x/E 

nential  particle  spectrum  at  high  energies  that  behaves  as  N E e 
On  the  basis  of  observational  indications  of  this  behavior  we  have  assumed 
6 = 2 as  a likely  value.  For  this  case,  we  may  also  rewrite  the  results 


(52)  and  (54)  in  a simple  form  in  terms  of  the  low  frequency  cutoff 
and  we  offer  the  following  speculations  on  second  stage  acceleration  in 


conclusion: 


The  choice  of  scaling  the  diffusion  coefficients  in  terms  of  the  quantity 
Emp  (21)  throughout  this  paper  was  motivated  by  the  possibility  of  deter- 
mining the  correlation  time,  from  observations  of  the  highest 
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1 


recorded  proton  energies.  The  value  taken  of  = lOGev  corresponds 


to  p = 4.3  Gev  and  from  (21)  we  infer  a correlation  time  and 


distance  of 


corr 


f lOOG  ]2 

f 1** 

^ ^mp 

^ 6x10^ 

lOGev 

sec 


X = 1.5 
corr 


f lOOG 

r E \ 

mp 

lOGev 

(56) 


km 


It  is  possible,  however,  that  the  correlation  time  is  much  longer  than 
(56)  and  the  much  higher  proton  energies  possible  are  not  recorded  because 
of  the  rapid  decrease  in  the  spectrum  (40)  at  very  high  energies.  Treating 


uJq  then  as  a free  parameter  and  taking  n = I5,  we  may  rewrite  in  summary 


the  results  for  spatial  diffusion  as  [cg  = 

^ ~ 


nn-  8 “0  "B 


(57) 


and 


" -571^ 


o^B 


(58) 


1 


Defining  a mean-free-path  such  that  K||  = y v x we  have  for  spatial 
diffusion 


o^B 


(59) 


The  energy  diffusion  coefficient  (19)  becomes 

D(n=^)  = J Ej  eJ 

leading  to  the  temporal  break  energy  (36)  of 

Eg(n=i5,t)  = ( ^ ] ^ ^^^■^corr^' 


(60) 


(61) 
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and  a steady-state  break  energy  for  diffusive  losses  of 

If  we  estimate  from  (16)  an  energy  diffusion  rate 

we  get  a result  which  differs  from  (61)  by  a factor  of  4 but  is  equal  to 
Melrose's  [1974]  result  [see  his  equation  24  and  the  discussion  following]. 
Equation  (62)  may  be  scaled  as 

Ep(n=i5)  = 18  f Y f f — ^ 1^  Mev/nucl  , (64) 


r 0.05s  Y 

■ 1 ] 

^ ^corr  ' 

1 1 J 

' lO^cm  - 

and  it  is  clear  that  longer  correlation  times  than  (56)  require  a larger 
Eg  in  proportion  to  achieve  a break  energy  in  the  range  of  20Mev/nucl. 

In  fact  if  we  consider  that  the  value  (64)  for  the  break  energy  is  typical 
we  may  place  the  correlation  time  in  the  range 

SxlO’-s  < < 5xl0'2s  , (65) 

where  the  smaller  limit  is  determined  by  the  parameters  given  by  (56)  such 
„ > 4.36ev  and  the  larger  limit  is  determined  by  (64)  such  that 

inaX  fP 

Eg  < 1. 

The  limit  eg  = 1 offers  an  interesting  interpretation  in  terms  of  low 
frequency  MHD  standing  waves  of  the  trap.  If  the  diffusive  scale  length 
is  L/2  and  the  wavenumbers  of  the  magnetic  perturbations  take  the  harmonic 
values  k = where  < L,  the  break  energy  for  diffusive  losses  (62) 

can  be  written  as 

"I  ( "corr 

where  k = 2irn  is  the  dominant  wavenumber  excited.  If  we  require  an 
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► 
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i 


1 

! 


energization  of  about  Eg  = 10^  E^,  and  suppose  that,  say,  L = 3X,  then 
the  limit  eg  = 1 is  consistent  with  a low  harmonic,  n 


corr 


6,  being 


excited  during  the  acceleration.  The  fundamental  frequency,  = 
2wV^/jJ,  gives  a greatest  upper  bound  to  the  correlation  time  of 


T < T 

corr  c 


r looG  1 

% 

1’’  f * 1 

1 K J 

■ 6xl0^cm"3 

■’  lO^cm  - 

(67) 


We  have  in  this  paper  considered  the  possibility  that  second  stage 
acceleration  may  occur  during  an/or  immediately  following  the  impulsive 
phase  of  proton  flares.  If  primary  energy  release  is  associated  with 
large  amplitude  Alfven  waves  self-consistently  generated,  then  (66)  admits 
the  possibility  of  low  harmonic  standing  waves  playing  a role  in  the  basic 
flare  instability.  Phenomenologically,  we  may  speculate  that  the  correla- 
tion time  (65)  is  related  to  the  pulse  width  of  hard  x-ray  bursts. 
Presently,  observations  show  a fine  structure  to  the  bursts  consisting 
of  -v  Isec  spikes  [Svestka,  1976].  The  Goddard  x-ray  spectrometer  aboard 
the  upcoming  Solar  Maximum  Mission  (SMM)  satellite  will  have  a temporal 
resolution  lO”^  s < t < 10"i  s from  20-300  Kev  energies.  If,  indeed, 
there  is  a relation  between  the  Alfven  wave  correlation  time  and  the  hard 
x-ray  pulse  widths,  we  believe  that  the  chances  are  good  that  this  instru- 
ment will  detect  fine  structure  in  the  range  of  (65).  The  University  of 
New  Hampshire/Max  Planck  Institute  gamma  ray  experiment  also  aboard  SMM, 
with  a temporal  resolution  down  to  64  ms,  should  settle  the  question  of 
how  soon  are  > 10  Mev  protons  accelerated. 
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FIGURE  CAPTIONS 


Figure  1.  The  temporal  evolution  of  the  hard  x-ray  and  2.2  Mev  v-ray 
bursts  for  the  August  4,  1972  flare  (from  Lin  and  Hudson, 
1976). 

I Figure  2.  Resonant  interaction  of  super-Alfvenic  protons  with  a linear 

polarized  elemental  scattering  center  of  length  x < 

The  solid  line  portions  of  the  trajectories  indicate  posi- 
tions on  this  side  of  the  page.  Energization  or  deenergiza- 
tion depends  on  the  cyclotron  phase  of  the  test  particle. 

An  average  over  a random-phased,  isotropic  distribution  of 
particles  leads  to  cyclotron  damping  of  the  oscillation. 
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